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Abstract 

We study the decay width of the inclusive process b — > sgg in the two Higgs dou- 
blet model with three level flavor changing neutral currents (model III). We analyse the 
dependencies of the differential decay width to the s- quark energy Es and model III 
parameters, charged Higgs mass mf^± and Yukawa coupling Cn bb- observe that there 
exist a considerable enhancement in the decay width for the relevant process. This en- 
hancement can be reduced by choosing C^^^ as negative and increasing the lower bound 
of mjj± to the large values, such as 800 GeV . This is an interesting result which gives an 
idea on the mass mfj± and sign of Cj^^ . 
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1 Introduction 



The B-meson system is interesting and rich phenomenologically, providing a comprehensive 
information on the theoretical models. With the forthcoming experiments at SLAC, KEK B- 
factories, HERA-B and possible future accelerators, the large number of events can take place 
and various branching ratios of events, CP- violating asymmetries, polarization effects, etc... 
can be measured ^. This will lead to test the models underconsideration and to determine 
corresponding free parameters. 

Loop induced processes are sensitive to the parameters of the models used. Therefore, they 
open a window for the determination of these parameters and investigation of new models. 
Among these type of decays, inclusive b —* sg reached a great interest since it is theoretically 
clean and sensitive to new physics beyond SM, like two Higgs doublet model (2HDM)|]^, mini- 
mal supersymmetric Standard model (MSSM) 0, |^, etc... . The Branching ratio Br oih sg 
decay in the SM is Br{h sg) ~ 0.2% for on-shell gluon and the enhancement of this ratio 
brings an advantage to decrease the averaged cham multiplicity rjc and to increase kaon 
yields 0. This enhancement can be obtained by including the QCD corrections or looking for 
new models beyond the SM. In the literature, there are number of theoretical calculations on 
the Br of the corresponding process beyond the SM. In |T^, |TT|, Br {b —>■ sg) was calculated in 
the 2HDM (Model I and II) for m/^± ~ 200 GeV and tan /9 ~ 5 and it was found that there was 
an enhancement less than one order. This decay was studied in the supersymmetric models 



1^ and further, the Br was calculated in the framework of the model III |jT3[, resulting with 
the enhancement at least one order compared to the SM one. This make it possible to describe 
the results coming from experiments fl^. 

In the case of time-like gluon, namely b —>■ sg* decay, Br should be consistent with the 
CLEO data [|ll 

Br {b ^ sg*) < 6.8% (1) 



and in |T^, it was showed that the model III enhancement was not contradict with this data 
for light-like gluon case. 

As a further process, g* can decay into quark-quark qq or gluon-gluon (gg) pairs. Inclusive 
three body decay b — > sgg is another interesting one which is studied in the literature extensively 
\[7\ , |TB|. It becomes not only from the chain process b —>■ sg* followed by g* — > gg but also 



from the emission of on-shell gluons from the quark lines to obey gauge invariance. In [|17 
the complete calculation was done in the SM and Br ratio was found at the order of 10~" 
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In |13| , |T8[ the additional contribution of gluon penguins in the Model III was estimated as 
negligible. 

This work is devoted to the study of the complete calculation for b sgg decay in the 
model III and we find that the decay width (F) is strongly sensitive to the charged Higgs mass 
mii±. Therefore, it is possible to get a considerable enhancement in F even 2 orders larger 
compared to the SM case. 

The paper is organized as follows: In Section 2, we give a brief summary for the model 
III. Further, we calculate the matrix element and decay width of the inclusive h sgg decay 
in the framework of the model III. Section 3 is devoted to discussion and our conclusions. In 
Appendix, we present the form factors appearing in the SM. 



2 The inclusive process b sgg in the framework of the 
model III 

The Yukawa interaction in the model III can be defined as 



= v7jQiL(l)lUjR + r]FjQiL4>lDjR + i^jQiL(t)2UjR + ^PjQiL4>2DjR + h.C 



(2) 



where L and R denote chiral projections L{R) = 1/2(1 =i= 75), 0j for i = 1,2, are the two scalar 
doublets. The Yukawa matrices rjfj^ and have in general complex entries. With the choice 
of 01 and 02, 



1 
72 







+ 



J_( V2H+ 



(3) 



and the vacuum expectation values, 

< 01 > = 



-r 



< 



>= 



(4) 



the SM particles are collected in the first doublet and particles due to new physics in the 
second one. The part of Yukawa interaction which is responsible for physics beyond the SM is 
the Flavor Changing (FC) interaction and can be written as 



where the couplings ^^'^ for the FC charged interactions are 



(5) 



VcKM 



(6) 



and ^^'^ is defined by the expression (more details see |]T9|) 

= (If . (7) 

Note that the index "N" in ^^''^ denotes the word "neutral". 
Now we start with the decay amplitude of the decay b sgg 

M{b - sgg) = ^^e^^{h)el{h)-s{p'K.'b{p) , (8) 

where e^(A;) are polarization vectors of the gluons with color a and momentum k. Using the 



rpah _ rp uT^ 



same parametrization for T^^; as in |T^, we have 

2 2 ''''2 2 

and can be obtained by the replacements ki ^-^ k2 and /i ^ z/ in the function T^^. Here ^ 
are the Gell-Mann matrices. The functions T^^^ and T^, in general, contain masses of internal 
quarks u,c,t in the SM and also d,s,b in the model III, since the process underconsideration 
takes place at least at one loop level. Therefore, at this stage, we take into account two different 
possibihties, 

• the mass of internal quark is heavy (namely, t-quark), 

• the mass of internal quark is light (namely, d, s, b, u, c-quarks). 

In the heavy internal quark case, the terms fcextemai/"^? ^external/ ''^1 (''^w^ ''^'h±) 
neglected and the form factors are obtained as functions of Xt = m^/m^r and yt = mf/m'jj 
where mfj± is the mass of charged Higgs boson in the model III. Neglecting s-quark mass, T^^ 
for the heavy internal quark is given by 

rpheavy _ • \ jp2HDM ^ / '^P'v'^ If 1^2 ^ La , ^ ^I'-ffj, ^ 

+ ^{2aai3kik§gf,u + '2ai,ak2^,q" -2af,akiuq" + (rf,uq^)}mbR . (10) 

Here q is the momentum transfer, q = ki + k2, Xt is the CKM matrix combination Xt = VtbV^* 
and p^H^^ is the form factor 

p2HDM ^ pSM (^^^ ^ pBeyond ^^^^ ^^^^ 

where -^2^^^ {xt) is the magnetic dipole form factor of 6 ^ sg* vertex (see Appendix). p^'^y°"-'^ (yj 
is the contribution coming from the charged Higgs boson in the model HI: 



^ts ^tb 

1 T/* T/ 

J- /r-rr r^rr ''r.c\ /rn rn tf 



t Vts Vtb 

i^N,U + eSc^) + ilsb^-^)G2{yt) , (12) 

rritnib Vts nb 



and 



G2{yt) = ^ {yt (3 - + yl) + 2 (-1 + y,) y, Iny,) . (13) 



In eq. (|T2[) we used the redefinition 



e--(-^i"-. (14) 

Note that we neglect the chiral partner of the form factor p^^y°^'^ (jj^^ and the neutral Higgs 
boson effects which should be very small due to the discussion given in pO| (see also Discussion 
part). 

If the internal quark is light {u or c) , the first additional contribution comes from mf / m^r 
and mj/mj^± terms. In the approximation 'rnf / 0, it is enough to replace ^2^*^ (xj 
with " — -F2(0)" since Ac = —Xt by unitarity, namely J2i=u,c,tXi = 0- There is no additional 
term coming from light quark for -^2^^^°"''^ (?/t) , since i^^^^?'""'^ (Q) almost vanishes. For light 
internal quark, the second additional contribution comes from k1xternail''^'i term which can 
not be neglected as in the heavy internal quark case. This contribution was calculated in the 



literature 17] and we give its explicit form in Appendix. Therefore, the resulting amplitude 



can be written as 

rp rpheavy i rplight . rplight M 

and T^^°'^y is given in eq. (0), T^i*^^* can be obtained from T^^'^'"^ with the replacement 
p2HDM _ -Fi^{Q) and T^jl' is given in Appendix. 

The decay amplitude for the process h sgg, Tl^^, can be parametrized by seperating color 
symmetric and antisymmetric parts [|l^ as 

-A^ A", _ ,A^ A", 

1 + 

2 ' 2 ^ ^.^[ 2 ' 2 



V' = TU^, -} + T-,[-, -] , (16) 



with 



1 

1 , 



= ^{T,u-Tp. (17) 
Finally we get the differential decay width of the process using the expression 



where Eg is the s-quark energy and Ei is the energy of photon with polarization e^(A;i). Here 
M is the average decay amphtude M = 2jip[ w ^ J = ^, Nc = 3. Now, we divide the 
differential decay width into sectors as follows: 

• Symmetric sector, (F"^^™), 

• Antisymmetric sector, (j'^^v"^^^ 

or 

• Right sector, (T^), 

• Left sector, (F^), 

• Left-rigth mixed sector, (F^^). 

Antisymmetric and symmetric sectors do not mix and they enter into decay width as 

j.Sym (Asym) _ Tr{T+'^-^ + mfe)) f /) P^^^' P""' , (19) 

with the corresponding color factors C+ = 'l^j^" and C_ = ^'''■^'^ respectively. Here 
we choose the polarization sum of the on-shell gluons as 

pMM _ —g^^^^ _|_ - ^ ^ ^ ^ 



ki.k2 

and T'^il7^ = 7o V 7o- Right sector contains form factors which are functions of Xi = 
mf/m^ and yi = mf/rn^j where i = u,c,t. Left one have the form factors which are cre- 
ated by the nonvanishing kl^ternai/''^'iight terms. Left-right sector contains mixed terms and its 
contribution is negligible. 

Since there are collinear divergences at the boundary of the kinematical region, we follow 
the procedure given in Jl^, namely taking a cutoff c in the integration over phase space: 



-2Esmb + ml + ml -2 E, mi, + ml 

+ cmt, < El < — - — — , (20) 



2 mfe 2{2Es-mb) 

with c = 0.01. Note that these limits are used in the integration over Ei and to get differential 
decay width j^. 
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3 Discussion 



There are many free parameters in the model III such as Yukawa couphngs, where i, j 
are flavor indices, masses of charged and neutral Higgs bosons. The procedure is to restrict 
these parameters using the experimental measurements. Since the contributions of the neutral 
Higgs bosons Hq and Aq to the Wilson coefficient C^-^'^ should not contradict with the CLEO 
measurement 



Br{B X,7) = (3.15 ± 0.35 ± 0.32) 10" 



(21) 



the couphngs ^§is{i = d, s, h) and S,^ should be negligible (see for details). In addition, the 
constraints |22|, coming from the AF = 2 mixing, the p parameter |2^, and the measurement 
by CLEO Collaboration results in the following restrictions: ^N^tc « ^nuj ^Nbb ^^^^ ~ 
0,C,^ij ~ 0, where the indices i,j denote d and s quarks. Therefore, we can neglect all the 
couplings except ^^^^ and ^^^b- This leads to the cancellation of the contributions coming from 
primed coefficient p^^^y"'^ and from the neutral Higgs bosons ho and Aq, having interactions 
which include the Yukawa vertices with the combinations of and ^§ ^g. Finally, we only take 
into account the multiplication of Yukawa couplings, ^^,tt^N^bb \^N,tt\'^ expressions. 

In this section, we study the s quark energy Eg, Yukawa coupling ^^^f, and charged Higgs 
mass mH± dependencies of the differential decay width for the inclusive decay b sgg. 
In our analysis, we restrict the parameters ^^u^ ^mb using the constraint for \Cy^-^\, 0.257 < 



IC^^^I < 0.439 [|T1, where the upper and lower limits were calculated in following the 



procedure given in Throughout these calculations, we take the charged Higgs mass m/^± 
400 GeV, and we use the input values given in Table dll). 



Parameter 


Value 


rric 


1.4 (GeV) 


rUb 


4.8 (GeV) 


Xt 


0.04 


rrit 


175 (GeV) 


mw 


80.26 (GeV) 


mz 


91.19 (GeV) 


J^QCD 


0.214 (GeV) 


as{mz) 


0.117 


c 


0.01 



Table 1: The values of the input parameters used in the numerical calculations. 
In Fig. 0we plot with respect to the s quark energy Eg, for ^^^b = 40 m;,, and \rtb\ 
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\-nr^\ < 1- tE~ is restricted in the region bounded by dotted (dashed) hnes for C^'^^ > 
{C^"^-^ < 0). Sohd hne represents the SM contribution. There is a considerable enhancement in 
the differential decay width especially for C^^^^ > case. Besides, the allowed region becomes 
larger for C^^-^ > 0. 

Fig. |]is devoted to the Eg dependence of color antisymmetric and symmetric part of j^. 
The color antisymmetric part lies in the region bounded by dash-dotted (dotted) lines and 
the color symmetric part by dashed (solid) lines, for C^^^ > {Cj^^ < 0). This figure shows 
that, for C^^^ > 0, the contribution of the color antisymmetric part is greater than that of 
color symmetric one. This is true also for Cj-^'^ < case. However the contribution of the 
color symmetric part for Cj^-^ > exceeds that of the color antisymmetric one for C^^-^ < 0. 
The allowed region becomes narrower for 0^"^-^ < (see dotted and solid lines). Note that the 
contributions due to the SM is presented by the solid and dashed lines which almost coincide 
with the X-axis. 

Fig. 1^ shows the Eg dependence of right, left and left-right mixed parts of -j^ in the SM. 
Solid line represents right, dashed line left and dotted line left-right contributions. The left 
one exceeds the right one up to almost Eg = 1.6 GeV since the kl^temai/^Ught contribution, 
responsible for left part, is comparable with the heavy internal quark, namely rrit, contribution. 
Left-right mixed part is very small and has also negative values. For the model III we have only 
right additional contributions since there is no left part beyond the SM in our approximation. 

In Fig. ^ and Fig. |^, we present the ^^^^ dependence of of for fixed values of -Ei = 2 GeV 
and Eg = 1 GeV. It is seen that there is almost no dependence on the parameter ^§ f,i, especially 
for its large values. 

For completeness, we also present rriHi dependence of for fixed values of ^^^^ = 40 m;,. 
El = 2 GeV and Eg = I GeV for G^^^ < (Fig. |). Here the restricted region is bounded by 
solid lines. This figure shows that there is a strong dependence on the mass m/^±. 

Now we would like to give some numerical results for our calculations. The total decay 
width for h sX transition is 

Ttot = (rlKbP +s|K,|2)ro , (22) 

where Fq = -^|^r and r, s are QCD sensitive parameters 

6.46 < r < 7.55 , 
2.38 < r < 2.92 

for ag = 0.2. Eeasy calculation shows that the total decay width is Ttot = 3.50 ± 1.50 10^^^. 
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In our calculation, we obtain the decay width for the SM as Tsm = 5.1 lO^^^GeV. For 
mH± = 400 GeV and ^^ ^^ = 40 rrib, the model III result is four (three) orders larger for C^'^-^ > 
{C^-^^ < 0) compared to the SM result. This is a strong enhancement contradict with the total 
decay width given above. This forces us to choose the sign of C^-^^ as negative {C^-^^ < 0) and 
also to take large values of charged Higgs mass, uih^ ■ For mn^ — 800 GeV, ^^^^ = 40 rub and 
Gj^^ < we get: 

8.58 10-^^ GeV < F < 1.5 10"^^ GeV , 
3.09 10"^^ Gey < F^2/rn < 2.70 lO^^^Gey , 

5.49 10"^^ GeV < T^^^"' < 1.23 lO'^^ GeV , (23) 
8.0810-^^^61/ < F^ < 1.45 10-^2 Gey , 
4.93 10-^^ GeV < F^ < 4.93 10"^^ GeV . 

In conclusion, we get a considerable enhancement in the decay width of the process b — > sgg 
in the model III. The enhancement can be suppressed by choosing G^^-^ < and increasing 
lower bound of charged Higgs mass, niHi- Further, the decay width of the process under 
consideration is not sensitive to the parameter ^^^b- 
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Appendix 



A The form factors in the SM for b sg* decay 

Here we present the magnetic dipole form factor Fg^*^ (xt) and the additional form factors due to 
the non- vanishing k'^xtemai/^Ught terms, (for details see [jl^])- The vertex function for b — > sg* 
decay with on-shell quarks can be written as 

^^l{p, P', q) = Fi {xt) (g%^ - 4) L - F2 {xt) i q" {nib R + rrisL) , (24) 

where p, p' and q are four-momentum of 6-quark, s-quark and gluon respectively. The magnetic 
dipole form factor ^2^^^ [xt) in the SM is 

pSMf^ ^ _ -8 + 38 - 39 xf + Uxf -5xj + 18 xf In xt 

12(-l + xt) 



and Xt = mf/m^. The non- vanishing k'^xternai/''^'iight terms for light quarks bring new additional 
contributions, AFi, Ai2, and A is (See Jl^ for details): 



where 



9 3 z 3 

9 z 3 z 3^^ ' ' 

At, = -1-2^^ , (26) 

z 



11— 

Qo(z) = -2 - iu+ - uJ)iln hivr), 

u+ 

1 ?/ 

Q.{z) = (Zn— + ^v^)^ (27) 

2 Ua. 



with 



w± = i(l± Jl-^) , (28) 



and 



g2 

z = — ^, i = u,c . (29) 
Finally, the contributions due to the non- vanishing k'^xtemai/'^'iight terms are 

- 2AFi(7,A;2^-7m^i.)^} (30) 
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Figure 1: a function of Es for fixed = 40 and |rt6| = Ith^I < 1- Here is 

^ ^N,bb ^ 

restricted in the region bounded by dotted (dashed) hues for C^^^ > {€7^^ < 0). Sohd hne 
represents the SM contribution. 
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Es {GeV) 



Figure 2: The color antisymmetric and symmetric part of as a function of Eg for fixed 

^Nbb ~ 40mfe, mH± = 400 GeV and |r(fe| < 1. Here the color antisymmetric part lies in the 
region bounded by dash-dotted (dotted) lines and the color symmetric part by dashed (solid) 
lines, for Cj^^ > {C^^^ < 0). The SM contribution almost coincides with x-axis. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

(GeV) 



Figure 3: Right, left and left-right mixed parts of a function of Eg for fixed ^^^^ = 40 m^, 
mH± = 400 GeV and \rtb\ < 1- Here sohd line represents right, dashed line left and dotted hne 
left-right contributions. 
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Figure 4: 

Cf^ < 0, El 
Dashed line represents the SM contribution 
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Figure 5: The same as Fig. ^ but for Cj^^ > 0. Dashed hne represents the contribution for 

C'," < 
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